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TECHNICAL NOTE 3690 

NORMAL COMPONENT OF INDUCED VELOCITY IN THE VICINITY 
OF A LIFTING ROTOR WITH A NONUNIFORM DISK LOADING 
By Harry H. Heyson and S. Katzoff 

SUMMARY 


Part I presents a method for computing, from available calculations 
for uniform disk loading, the effect of nonuniform circularly summetrical 
disk loading on the normal component of induced velocity in the vicinity 
of a lifting rotor. Charts of the normal component of induced velocity 
are given for the longitudinal plane of symmetry and for the major axes 
of rotors with two different, nonuniform circularly symmetrical disk 
load distributions. It is shown that the normal component of induced 
velocity must be zero at the center of any practical rotor. 

A comparison of the results of this paper with those for a uniform 
disk loading shows that nonuniform disk loading has a powerful effect 
on the induced velocity distribution and that it must be taken into 
account in estimating the effect of the rotor on most components of an 
aircraft. 

Part II develops certain symmetry relations for the induced veloc- 
ities in the plane of a uniformly loaded rotor and also develops relations 
between the radial load distribution of the rotor and the radial varia- 
tion of induced velocities in the wake. 


INTRODUCTION 


Recent rotary-wing designs incorporate tails for increased s+-.bility 
and auxiliary wings for improved forward-flight performance. Any estimate 
of the behavior of the complete aircraft depends upon a knowledge of the 
flow induced by the rotor in the neighborhood of these auxiliary lifting 
surfaces. Such information is relatively meager. The only available 
analytical treatment is that of reference 1 , which calculates the normal 
component of induced velocity of a uniformly loaded lifting rotor along 
its major axes and in its longitudinal plane of symmetry. Although the 
assumption of uniform disk loading and the fact that the calculations 
were made for only one plane are obvious limitations of reference 1 , it 
was nevertheless hoped that the results would give useful indications 
of the downwash over small-span auxiliary lifting surfaces. 
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Wind-tunnel flow surveys (ref. 2), however, indicate that there may 
be very little relation between the downwash given in reference 1 and the 
downwash for an actual rotor. For example, large variations in induced 
velocity across the lateral axis of the rotor were observed, whereas the 
calculations show no variation at all across the lateral axis. Probably 
the main cause of the variation of induced velocity along the lateral 
axis is that the disk loading is not uniform as is assumed in reference 1. 
The purpose of the present paper is to study the consequences of this 
nonuniformity of disk loading. 

Part I describes a method of superposition, by means of which the 
results of reference 1 may be used to calculate the induced flow about 
rotors having an arbitrary circularly symmetrical distribution of rotor 
disk loading. Calculated results are presented for two rotors having 
different nonuniform disk loadings. These results show how the distribu- 
tion of disk loading affects the downwash, and are used in reference 2 to 
explain the measured distribution of induced velocity. 

Part II contains analytical treatments of some interesting charac- 
teristics of the downwash field that were observed during the course of 
the calculations for part I. More specifically, part II develops certain 
symmetry relations for the induced velocities in the plane of the uniformly 
loaded rotor and also develops relations between the radial load distribu- 
tion of the rotor and the radial variation of induced velocities in the 
wake . 


I. DEVELOPMENT OF INDUCED VELOCITY FIELDS BY SUPERPOSITION 

By Harry H. Hey son 

Symbols for Part I 

b number of blades 

L lift, lb 

r blade-element radius, ft 

R rotor radius, ft 

Ry radius of vortex cylinder, ft 

v normal component of induced velocity, positive downward, ft/sec 

Vq average, or momentum, value of the normal component of induced 

velocity, positive downward, ft/sec 
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X distance parallel to longitudinal rotor tip-path-plane axis, 

measured positive rearward from center of rotor (fig. l) , ft 

Y distance parallel to lateral rotor tip -path-plane axis, measured 

positive on advancing side of disk from center of rotor 
(fig. 1), ft 

Z distance parallel to vertical rotor tip-path-plane axis, 

measured positive upward from center of rotation (fig. l), ft 

T circulation, ft^/sec 

T* circulation at tip of blade, ft^/sec 

p mass density of air, slugs/ft3 

X rotor wake skew angle, angle between the Z-axis of the tip-path 

plane and the axis of the skewed wake, positive rearward from 
Z-axis (fig. l) , deg 

rotor angular velocity, radians/sec 


Theory 

The wake considered as a system of coaxial vortex cylinders . - For 
the present simplified theory, the forward speed of the rotor is assumed 
negligible relative to the rotational speed of the blades. Then, at any 
radial location, the local velocity of a blade element is simply Hr, so 
that the lift per unit radial distance can be written 


-3— = bpQrf (1) 


where T, in general, is a function of r. The local disk loading (or 
annulus loading) is 


t , bpfir T dr opHT (o s 

Local disk loading = — = — — \2) 

2«r dr 2 * 

In general, as illustrated by equations (l) and (2), the local disk 

loading varies as — times the blade loading. Thus, a uniform disk 

loading implies a triangular blade loading, and a triangular disk loading 
implies a parabolic blade loading. 
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Equation (2) also shows that for a uniform disk loading the circula- 
tion is constant along the blade. For this case, vorticity is accordingly 
shed only at the blade tips, and the wake consists of a single vortex 
cylinder. In references 1 and 3 , where the lifting rotor is represented 
by a uniformly loaded disk, the downwash field is correspondingly calcu- 
lated as the downwash field of this single vortex cylinder. 

If, however, the disk loading is not uniform, the circulation varies 
along the length of the blade and vorticity must be shed all along the 
blade instead of only at the tip. Consequently, a rotor having nonuniform 
disk loading must be represented by a distribution of shed vortex 
cylinders . 

In the present calculations, these cylinders are all assumed parallel 
and concentric (see fig. 2(b)) and axe inclined at a skew angle defined 
by the forward velocity of the rotor and the momentum- theory value of 
the rotor induced velocity. Such an assumption presents obvious incon- 
sistencies, because such an idealized geometry could not be produced or 
maintained in the nonuniform induced-flow field that it itself creates. 

In spite of such inconsistencies, however, the arrangement of concentric 
cylinders is probably a better basis for computations than a single vortex 
cylinder trailing from the blade tips. 

It is possible to consider that the solution of reference 1 represents 
the flow generated by any one of these vortex cylinders. The radius used 
in that report should now be considered to be the radius R v of the 
individual vortex cylinder. The contribution of any single vortex 
cylinder to the total induced flow at a point ma y "be read directly 

from the charts of reference 1 at the point [(x/R) (R/R v ) , (y/r)(r/R v ), 
(Z/R)(R/R V )] . 

In practice, any arbitrary circularly symmetrical rotor disk- load 
distribution may be approximated by the sum of the loads carried by a 
finite number of superposed concentric vortex cylinders of suitable 
strength and dimensions. Then the induced velocity field of the non- 
uniformly loaded rotor will be the sum of the superposed induced velocity 
fields of the vortex cylinders comprising its wake. 


Some difficulty may be anticipated because of the discontinuities in 
the flow at the edges of the vortex cylinders. In practice, however, only 
a small error occurs if a sufficient number of vortex cylinders are 
used, and the points of discontinuity themselves are avoided. 


Sample calculation .- In order to demonstrate the method of super- 

( X Y Z ^ 

— = 0.2, — = 0, = 0.1 

R R R ' 

will be calculated for a rotor with a triangular disk loading and 
operating at a skew angle of 90°. 
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As shown previously, the local disk loading is proportional to T; 
therefore, for uniform disk loading the circulation is constant along 
the blade. If the disk loading is triangular, 


r = r* 


r 

I 


so that the local disk loading is 


Local disk loading = 7 

dn R 

and the average disk loading for the entire rotor is 

1 f R bpfi™ r„ 

Average disk loading = — 7 / -7 — R p2mr < ^ r 

i<Rr n 

_ bpar' 

3* 

Thus, for this case of the triangular disk loading. 

Local disk loading = -- X (Average disk loading) 

2 R 


Thus, the vortex field of a rotor with a triangular loading may be 
represented by one positive vortex cylinder with the same radius as the 
blade tips and a strength 1.5 times as great as a uniformly loaded rotor, 
plus ten equal negative vortex cylinders of smaller radii and a strength 
0.15 times as great as a uniformly loaded rotor. The sum of the resulting 
disk loads is shown in the following sketch and adequately represents a 
triangularly loaded rotor: 
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The contribution of each vortex cylinder may be found from fig- 
ure 4 (g) of reference 1. 


The computation may be carried out in tabular form as follows: 


(R v /R) 

(X/R) (R/R v ) 

(y/r) (r/r v ) 

(Z/R) (R/R v ) 

( v / v o> v 

1.00 

0.200 

0 

0.100 

1 . 15 ( 1 . 50 ) = 1.725 

• 95 

.211 

0 

.106 

1 . 15 (-. 15 ) = -.173 

.85 

.236 

0 

.118 

1 . 17 (-. 15 ) = -.176 

• 75 

.267 

0 

.134 

1 . 19 (-. 15 ) = -.179 

-65 

.308 

0 

.154 

1.20(-.15) = -.180 

.55 

.364 

0 

.182 

1.21(-.15) = -.182 

.45 

.445 

0 

.223 

l. 23 (-.l 5 ) = -.185 

.35 

• 571 

0 

.286 

1 . 24 (-.l 5 ) = -.186 

• 25 

.800 

0 

.400 

1.25(-.15) = -.188 

.15 

1.330 

0 

.667 

. 95 (-.l 5 ) = -.143 

.05 

4.000 

0 

2.000 

.26(-.i5) = -.039 




v/v c 

1 = ,^( v / v o) v = °*° 9 ^ 


It will be noted, in general, that the charts of reference 1 must 
be extrapolated to obtain the contribution of some of the smaller 
cylinders. Also, if the maximum loading does not occur at the tips, 
more than one positive vortex cylinder will be required to adequately 
represent the flow. 

Induced velocity at center of rotor .- Consider a rotor with a disk 
loading that is uniform everywhere except in a central "cut-out" portion 
where it carries no load. See the following sketch. 



The vortex field of this rotor may be represented by just two vortex 
cylinders of equal, but opposite, vorticity. Since the smaller cylinder 
has the same downwash field as the larger cylinder (except for dimensions), 
it follows that the induced velocity is zero at the center of the rotor ’ 

disk. The same result, of course, will follow whenever any rotor has 
zero disk load at its center. Since all practical rotors have some cut- 
out near their centers (due to hub, etc.), they all must have zero induced 
velocity at their centers. 
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This result is a special case of a more general theorem which is 
occasionally useful in comparing the results of theoretical calculations 
made for several different load distributions. Consider a rotor with a 
step loading composed of a number of uniform loads of different radii 
such as in the following sketch: 



The vortex field of this rotor may be represented by several vortex 
cylinders of different strength. A uniformly loaded rotor produces an 
induced velocity of Vq at its center. Therefore, the induced velocity 
at the center of the step-loaded rotor will be vq 1 (due to uniform 

load 1) plus vq 2 (due to uniform load 2) and so on. The final result 

is that the induced velocity at the center of the step-loaded rotor will 
then be the same as if the local load at the center existed over the 
entire disk rather than just at the center. 


Calculated Results 

Longitudinal plane of symmetry .- A number of nondimens ional charts 
of the induced velocity have been computed. These charts not only illus- 
trate the effect of a nonuniform disk loading, but they may also be used as 
working charts for the determination of the induced velocity at points 
in the field of practical rotors. The calculations were made for rotors 
having two different, circularly symmetrical, nonuniform disk loadings. 
These assumed disk loadings are shown in figure 3. The triangular 
loading was chosen since it is the simplest approximation to the actual 
rotor disk loading in all flight conditions. The variation of disk 
load labeled "typical load" was measured at a typical cruising condition 
during unpublished results from tests in the Langley full-scale tunnel of 
a rotor equipped to measure the rotor-blade pressure distribution. 

Figure 4 presents contour charts of the normal component of induced 
velocity in the longitudinal plane of symmetry for the rotor with a 
triangular disk loading for six different skew angles in the range from 
0° to 90°. Figure 5 presents similar charts for a rotor having the 
typical measured disk load distribution for three different skew angles 
between 63 . 43 ° and 90°, which bracket the skew angle at which this load 



8 


NACA TN 3690 


distribution was measured. Figures 4 and 5 may also be used to find the 
induced velocity distribution for skew angles between 90° and l 80 ° (the 
autorotation range) if the Z/R scale is multiplied by -1 and the chart 
for the supplement of the desired skew angle is read. 

Longitudinal axis .- Figures 6 and 7 present the induced velocity 
distribution along the X-axis (or longitudinal axis) of the rotors with 
the triangular and typical measured load distributions, respectively. 

lateral axis .- Figures 8 and 9 present the induced velocity distribu- 
tion along the Y-axis (or lateral axis) of the same two rotors. 

Figures 6 to 9 may be used for skew angles from 90° to l8o° merely 
by reading the curves for the supplement of the desired angle. 

Figure 8 indicates the order of accuracy of the calculations. Part II 
of this paper will show that the variation of induced velocity on the 
latera.1 axis must be triangular if the disk loading is triangular. 

Actually, figure 8 indicates a slight upward curvature. The difference 
between these curves and straight lines is, however, small and could 
probably have been avoided by using a somewhat larger number of cylinders 
to represent the wake. 

Other regions .- It would be interesting to carry out these calcula- 
tions over the lateral planes where auxiliary devices such as wings, tail 
surfaces, or propellers would be located on flight machines. However, 
such calculations would be dependent upon a knowledge of the flow field of 
the uniformly loaded rotor and this flow field has not yet been computed. 
The flow in these regions can be calculated, if required, from considera- 
tions in reference 1 and the present paper. 


Discussion of Calculated Results 

Longitudinal plane of symmetry .- Figure 10 presents a chart of the 
induced velocity in the longitudinal plane of symmetry of a uniformly 
loaded rotor at a skew angle of 63 . 43 ° (reproduced from ref. 1). 

The effect of nonuniform axisymmetric loading may be seen by comparing 
figure 10 with figures 4(d) and 5(a) . Since both of the assumed nonuniform 
loads are zero at the center of the rotor, the induced velocity at the 
center is zero in both cases, as opposed to a value of v/vq = 1.0 at the 

center of the uniformly loaded rotor. Both of the nonuniformly loaded 
rotors show an appreciable area of upwash just below and behind the 
center of rotation. They also show zero induced velocity further rearward 
on the center line of the wake. In this plane, the chart for the triangular 
load distribution shows a roughly triangular induced-velocity profile across 
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the wake, and the chart for the typical measured load distribution shows 
a more rounded profile across the wake with a larger area of zero induced 
velocity reflecting the effect of the cut-out region in the vicinity of 
the hub. In contrast, the results for the uniformly loaded rotor case 
shows the induced velocity building up continuously as the flow passes 
through the rotor and rearward in the wake until it reaches a uniform 
value of v/vq =2.0 near the trailing edge of the rotor. 

Longitudinal axis.- Figure 11 (obtained from figs. 6, 7, and 10) shows 
the calculated effect of rotor disk- load distribution on the induced- 
velocity distribution along the X-axis (longitudinal axis) of the rotor. 
This figure clearly shows the large effect caused by the requirement 
that the induced velocity be zero at the center of the rotor for each of 
the nonuniformly loaded cases. Each of the nonuniformly loaded rotors 
has an upwash just behind its center, whereas the uniformly loaded rotor 
has downwash values of v/v 0 > 1 in this region. The differences in. 
induced velocity between the nonuniformly loaded rotors are small along 
this axisj however, the rotor with the typical measured load distribution 
does have a somewhat greater upwash behind the center as a result of the 
larger cut-out region. 

* 

Lateral axis .- Figure 12 presents a similar comparison along the 
Y-axis (lateral axis) of the rotors at a skew angle of 90°* This figure 
shows that the lateral center-line distribution is closely similar in 
shape to the distribution of load on the disk. This similarity is 
discussed further in part II. 

It should be noted that the distribution given for the uniformly 
loaded rotor is not the same as that given in reference 1, since in the 
present paper the induced velocity is shown as being uniform across the 
lateral axis. A proof that the induced-velocity distribution must be 
uniform at this location is given in part II. In discussions between 
the senior author of reference 1 and the authors of the present paper, 
it was suggested that the inaccuracy in the plane of the disk resulted 
from a basic difficulty in the computational procedure in the location 
near the edge of the disk. It is assumed that the yalues given in 
reference 1 are correct for locations outside the rotor tips. The 
calculations made for the present report were based on composite curves 
similar to the one shown in figure 12. 

An interesting observation can be made from figure 12. Reference 4 
has shown that the uniformly loaded rotor values of induced velocity 
in the plane of symmetry (ref. 1) can be used as a guide to estimate the 
average induced velocity across a span of the order of the rotor diameter 
at any location. The downwash values for the lateral axis justify this 
result to the same order of accuracy as the measurements of reference 4, 
since the average induced velocity across the span for both nonuniform 
disk loadings is approximately the same as that for the uniformly loaded 
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rotor. (That is, the average value of v/v Q is about 1.0 in all three 
cases.) However, the differences between the results of reference 1 and 
this paper in the longitudinal plane of symmetry indicate the possibility 
of large errors if the .uniformly loaded rotor flow field is used to 
estimate the effect of the rotor on items of appreciably different span 
than the rotor. 


Concluding Remarks for Part I 

A method of calculating the effect of nonuniform circularly symmet- 
rical disk- load distributions on the normal component of the induced 
velocity of a lifting rotor has been presented. 

The induced velocity at the center of any rotor which has zero 
load at its center (including all practical rotors) must be zero. 

Charts have been presented for the normal component of induced 
velocity along the major axes and for the longitudinal plane of symmetry 
over a wide range of skew angles for rotors with two different, nonuniform 
axisymmetric disk loadings that are representative of the actual loading 
on the rotor disk. 

Rotor disk- load distribution has a large effect on the induced- 
velocity distribution and must be taken into account in esti ma ting the 
effect of the rotor on most components of an aircraft. 


II. SYMMETRY RELATIONS AND THE RELATION BETWEEN 
RADIAL LOAD DISTRIBUTION AND THE RADIAL 
DOWNWASH DISTRIBUTION IN THE WAKE 
By S. Katzoff 


Symmetry Relations Concerning the Induced Velocities 
in the Plane of the Rotor 

The first section of pant II is concerned with the induced field 
of the skewed-cylinder vortex that is assumed in references 1 and 3 to 
represent the wake of a uniformly loaded disk. It will be shown, without 
reference, to the detailed equations of the flow, that the induced- 
velocity field produced in the plane of the disk by this idealized wake 
possesses certain symmetries. By providing relations between the induced 
velocities at pairs of symmetrically located points, such symmetries are 
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useful for checking the accuracy of calculated values such as those of 
reference 1. They also provide directly the values of induced velocities 
along the lateral axis of the disk. 

It will he noted that the terminology followed in this part of the 
paper differs from that usually used in rotary-wing work. The present 
notation was used in order to retain the usual complex-variable symbols 
in two-dimensional flow analyses and also to avoid a confusing multi- 
plicity of subscripts. The reader should note in particular that the 
normal component of induced velocity w is positive upward in the 

following material. 

Points on the rotor disk .- In figure 13(a) consider a pair of 
points on the rotor disk, such as P and Q, that have the same lateral 
location and are symmetrically located with respect to the lateral axis. 
That is, if P is at point (x,y,0), then Q is at point (-x,y,0). Let 
the three components of the induced velocity at P be (u^, v^, w^) and 

at Q be (U2, Vg, Wg) , where the third component is generally the only 
component of interest. 

Now consider the entire diagram rotated l80° about the Y-axis so 
that it appears as in figure 13(h) . Points P and Q thereby move to 
points P' and Q', where the induced velocities are, respectively, 

(-u x , v x , -w-l) and (-u 2 , Vg, -Wg) . 

If figure 13(b), with the vorticity reversed, is fitted to fig- 
ure 13(a), so that P* falls on Q and Q' falls on P, the two semi- 
infinite vortex cylinders form a single continuous infinite vortex 
cylinder (fig. 13 (c)). The induced velocity at the left-hand point 
(P' and Q) is (u^+Ug, Vg-Vj, w^+wg) , and at the right-hand point 

(Q 1 and P) it is (u^ug, v^-Vg, w^+wg) . The two velocities are thus 
equal except for the lateral components, which are equal and opposite. 

In appendix A of reference 3, however, it was shown that the induced 
velocity within this infinite skewed helix is uniform and has no lateral 
component. The first result of the present discussion, then, is that 
Vg ” v i = Vi - vg = 0 or vj = vg. 

The second result is that the sums u-j_ + Ug and w^ + Wg are 

uniform. That is, the sums of the two longitudinal components and of the 
two vertical components are the same for all pairs of symmetrically 
located points on the disk (as P and Q) , and are equal to the longi- 
tudinal and vertical components of the induced velocity within the 
infinitely long helix. In the nomenclature of reference 3> 

U 1 + Ug = v x ' 
w ! t wg = v z ' 
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Furthermore, if points P and Q are on the lateral (Y) axis, so 
that they coincide (u 1 = U2 and w-j_ = w 2 ) , 



That is, the longitudinal and vertical components of the induced velocity 
are constant all along the lateral diameter of the disk and are equal 
to the values at the center of the disk. 

Points outside the rotor disk .- If points P and Q are beyond 
the edge of the rotor disk (fig. 14), the discussion proceeds as before 
except that it no longer follows that u-j_ + u 2 and w-|_ + w 2 are con- 
stant or that vj - v 2 is zero. The reason, of course, is that after 

superposition of the two semi-infinite skewed vortex cylinders the 
points lie outside the infinite elliptic cylinder instead of inside of 
it. The combined flow field outside the cylinder is not uniform but is 
the same as if the cylinder had a solid boundary. That is, the com- 
bined induced flow is two dimensional in planes normal to the cylinder 
axis and can be computed as that induced by an elliptic cylinder in a 
steady flow of velocity V cos X normal to its axis. The streamlines 
of the induced flow in the normal plane are sketched in figure 15. The 
components u' and v' of the induced flow in this plane are related 
to ui + U2, vi - V2, and wq + W2 as follows: 


U 1 + u 2 = v' cos X 
wi + w 2 = v' sin X 

v 2 - v 2 = u' at the right-hand point (P and Q') 
For points along the lateral axis., 

U 1 + U2 = 2 u i = v' cos X 
w x + w 2 = 2 w x = V* sin X 
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Thus, for points along the lateral axis, and w^ can be determined 
directly from the known two-dimensional flow about the ellipse. 

In order to facilitate application of these results, a few remarks 
concerning the calculation of v' are contained in the next section. 

Flow about an elliptic cylinder .- Consider (fig. 15) the flow 
induced by an elliptic cylinder of semimajor axis a and semiminor 
axis b in a cross flow of velocity V cos X. Let the plane of the 
figure be the z' plane, where z' = x' + iy' (the x' axis is parallel 
to the y-axis, or the lateral axis of the disk). 

The flow is probably best treated in elliptic coordinates, related 
to the rectangular coordinates by 

z ' = c cosh £ 


where 

5 - i + in 

c = \Je£ - b^ 


The complex flow function w for the field induced by the ellipse 
in the steady cross flow of velocity V cos X is 

w = (V cos X) ia ^ ~ e“^ 

(The equation may be obtained, for example, from the equation on page 256 
in ref. 5 by putting U equal to zero, replacing V by V cos X, and 
omitting the free -stream flow function.) Then the complex velocity is 


dw _ dw/d£ 
dz • ~ dz ' /dt, 


-(V cos 



a + b e~£ 
a - b c 'sirih 


i 


- - (v cos x) t\/Ht pttt 

Putting C = i + iT l 811(1 multiplying numerator and denominator by the 
conjugate of the denominator - l) leads to the following expres- 

sion for v* , the negative of the imaginary part of 
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- 2 < v c ° s x > f /i? 


*26 


cos 2n -1 


e 1 ^ - 2e 2 ^ cos 2^+1 


or 


v.' = (V cos X) a ( a t b ) 
" - c 2 


cos 2ri - e~ 2 ^ 
cosh 2£ - cos 2 t) 


Since x + iy c cosh (6 + iT]), the values of | and T) corresponding 
to any point (x*,y*) in the field of the ellipse are obtainable directly 
from tables of complex hyperbolic functions (ref. 6). (in using these 
tables, note that t) is given in units of n/ 2. It will be observed 

that cos X = — and sin X = — 


It will be noted that the points P, P», Q, and Q 1 He in a 
plane which intersects the plane of the ellipse at the angle X. However, 
because of the two-dimensional character of the flow, the desired values * 
of the induced velocities at a point (x,y) in the plane of the rotor may 
be found at the point (x'=y, y'=x cos x) in the plane of the e lli pse. 


Relations Between the Radial Variation of Disk Loading and the 
Radial Variation of Downwash Velocity Within the Downstream Wake 


It was shown in figure 12 that the computed lateral variation of the 
downwash angle across the lateral axis resembled the assumed radial dis- 
tribution of disk loading. An effort was therefore made to determine 
whether the resemblance was mainly fortuitous or whether a general 
theoretical basis existed for it. 

It was found possible to show that, for the far wake, the downwash 
velocity along any radius is proportional to r n if the disk loading 
is also proportional to r n . The proportionality factor, however, 
depends on n (in addition to the azimuth). Accordingly, if the radial 
disk loading is represented by a power series, the downwash velocity along 
any radius is also represented by a power series where, however, the 
respective coefficients of the two series are not proportional. Thus the 
radial loading is not, in general, exactly proportional to the radial 
downwash— velocity distribution, although the two may be nearly propor- 
tional if the radial loading is approximately proportional to, say, 
r or r 2 . 

In the following development the simplest case, X = 0, will be 
discussed first. The general case, 0° < x < 90°, will then be discussed. 
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Finally, the particular case X = 90°> which is related to linearized 
lifting-surface theory, will he discussed. 

The case X = 0° .- The problem is simple for the circular wake, 

X = 0^ Consider first the uniformly loaded disk, which produces a 
single helical vortex lying along a circular cylinder. In the far 
wake the downwash within the cylinder is uniform, and there is no down- 
wash at all outside the cylinder. For an infinitesimal annulus (dr) of 
disk loading, the wake is then two concentric circular cylinders having 
equal and opposite vorticity and separated by dr. The downwash between 
the two cylinders is uniform and proportional to the annulus loading. 

The downwash everywhere else is zero. If the disk loading is made up 
of a continuous distribution of such annuli, the downwash directly 
downstream of each annulus is thus determined only by that annulus loading 
and is not affected by the loading of any other annulus. For this case, 
then, the radial variation of downwash velocity is proportional to the 
radial variation of disk loading. Restricting the theorem to loadings 
that are proportional to r n is unnecessary for this case, which will 
be recognized as merely the well-known propeller blade-element theory 
applied to the far wake instead of to the disk itself where the induced 
velocities are only half as much. 

The case 0° < X < 90 0 .- As previously mentioned, the uniformly 
loaded disk produces an ultimate wake in the form of an inclined elliptic 
cylinder, within which the induced flow is uniform and outside of which 
the induced flow is that due to an elliptic cylinder in the cross-flow 
component of the free-stream flow. If the loading is circularly symmet- 
rical but not uniform, the wake consists of a corresponding distribution 
of concentric similar elliptic cylinders, and the present problem is 
concerned with the superposition of their induced flow fields. 

If each cylinder induced only the uniform internal velocity field 
and induced no external field, the argument would proceed just as for 
the case X = 0°, and it would follow immediately that the radial distri- 
bution of induced velocity is proportional to the radial loading distri- 
bution. In order to prove the theorem, then, it is only necessary to 
show that the total of the external fields of the distribution of cylinders 
will, along any radius, also have velocities proportional to r n , if the 
loading is proportional to r n . 

Consider the flow induced by an ellipse of given shape and unit 
size (say unit semi major axis) moving at velocity V cos x. As previously 
indicated, the induced velocity component parallel to the minor axis 
(parallel to the direction of motion of the ellipse) is the component 
of interest. Along any direction 0 from the origin, this induced 
velocity is a function of the radial distance p. 
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If, while the shape remains constant, the linear dimensions of the 
ellipse are changed by the factor k, the velocity originally existing 
at radial distance p will be found at radial distance kp along the 
given direction 9. For the general ellipse of size k, then, the 
velocity component v* can be expressed as 

v' = f(0,k/p) 


The velocity induced by an increment in size, or annulus, corresponding 
to dk is 


^ak-f'(e,k/p)f 


where f' denotes the derivative with respect to k/p. If each incre- 
mental velocity is weighted by the factor k n , the total velocity is 



(k/p) n f(e,k/p)^ 


where the upper limit K is that value of k for which the ellipse just 
touches the specified point ( p, 9) . . (As previously noted, the present 
discussion concerns only the contributions of those annuli for which the 
point is located externally.) Changing the variable of integration 
from k to k/p changes the form of the preceding expression to 

r K/p 

fP (k/p) n f»(e,k/p)d(k/p) = p n F(9,n) 

J k/p=0 

since the upper limit K/p is now a function only of 9. 
induced velocity along any radial line is proportional to 
was to be proved. The proportionality: factor is a function of both 9 
and n. It will be observed that this part of the proof can apply to 
any shape, since the fact that the wake is elliptical was not required 
or used in the proof. 

Performing the indicated integration in order to get the proportion- 
ality factor F(9,n) is obviously a somewhat awkward task. Possibly 
the most convenient method is first to integrate by parts in order to 
return f ' to f (that is, v', for which an expression was previously 
presented) , and then determine the resulting integral numerically. For 
either the lateral or longitudinal axis, however, the calculation may 
be performed without excessive difficulty. The method will be here 
indicated. 


Thus , the 
p n which 
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For the lateral axis (t| = 0 ) , the previously derived expression for 
v' reduces to 


v' 


= ( v cos x ^Vr? 


b xj 
b ~ 


- 

x/x ' 2 - o 2 


so that the problem is reduced to evaluating 



,n 




Substituting c = x' sin a changes the integral to a form that can be 
integrated by application of items 27^- and 263 of reference 7* 

For the longitudinal axis t) = ^ the expression for v' reduces to 


-(V cos X) 


'\£2 


+ y' d - y’l 


\/c 2 + y 2 


so that the problem is reduced to evaluating 


b y, c n aj \/c2 + y ,g - y 


C=0 


7 


dc 


y“ 


Substituting c = y 1 tan a changes the integral to a form that can be 
evaluated by application of item 27^ of reference r J. 


The case X = 90 ° .- The limiting case of the flat wake is frequently 
assumed for convenience in analyses pertaining to high-speed forward flight. 
Since the rotor and the wake are in the same plane the case is analogous to 
that of classical linearized wing theory, with the wing in this instance 
having a circular plan form and a circularly symmetrical loading. 


Consider first a uniformly loaded rotor. The "span load" distribu- 
tion of the disk, considered as a wing, is then elliptical, since it is 
merely proportional, at each spanwise station, to the "chord". As is 
well known, the downwash in the far wake behind an elliptically loaded 
wing is uniform. That is, disregarding the rolling-up phenomenon, the 
flow in the vertical plane normal to the far wake is the two-dimens.ional 
flow about a horizontal straight line of span equal to the rotor diameter. 


moving downward with a velocity equal to 


2 C l V 2 C l V C l V 

rt (Aspect ratio) 2 
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(see fig. l6). The complex velocity for this two-dimensional flow 
around the far wake is given by 


where 

w 



iv* 



i 


complex flow function 


r semispan of wing (or radius of the rotor) 

z 1 = x 1 + iy 1 where , in the present case (X = 90°); x * is the same as 
lateral dimension y and y 1 is the same as vertical 
dimension z 


If the radius of ‘the uniformly loaded circular wing is increased 
hy dr, the complex velocity in the flow field changes by 



This expression thus gives the complex velocity in the far wake contributed 
by a wing in the form of a circular annulus of radius r and width dr, 
and having a lift coefficient C^. 

Now consider a circular wing on which the loading is proportional 
to the nth power of the radius. It is made up of annuli of variable 
radius r and lift coefficient given by, say, A r n , where A is a 
constant. The total complex velocity in the far wake is found by 
substituting Ar n for Cl in the preceding expression and integrating 

with respect to r between 0 and the outermost annulus, r = R. The 
complex velocity is thus 


f R Ar n V 

-irz* 

4=o 2 

(z 1 ^ - r ^5/2 


dr 


where the upper limit R is greater than | z’| . 

In order to avoid a singular point in the integration, the point 
z f will be assumed to lie slightly above the real axis and will be 
allowed to approach the real axis after the integration is performed. 
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Substituting r 
the integral to 


z' sin 9 (where 9, in general, is complex) transforms 

R 


- 1 — 

-lAz«"V /' Sin " Z ' sin n+I e d9 
9=0 


1 ,n v r 

^ p_> n (i An2 


cos" 1 © 


vhich can be evaluated by first applying the second form of item 274 of 
reference 7 and then reducing jfsin n+1 9 d 9 by successive applications of 
item 263 . Except for the last term, the. integral is 

-iAz ,n vjsin n+1 9 tan 9 + sin n 9 cos 9 + sin n-2 9 cos 9 + 


n(n - 2 ) 

(n - l)(n - 3 ) 


sin n “^ 9 cos 9 + 


•••] 


Substituting the limits and letting z' become a pure real, less then R, 
results in a pure real; that is, all of these terms contribute nothing 
to the vertical-velocity component. The final term of the integral 
involves 


. -iL 

,sm z» 


sin 0 d 8 or 


• .i-R 

,sm Jm z , 


d9 


f 0 


'0 


depending on whether n is even or odd, respectively. For n even, 


R.P. 


1 R 


r- 

-i 

psin z ' 

4 

sin z ' 



/ sin 9 d9 = R.P.(- cos 9) 

> 0 

= R.P. 

o 

;\ 

i-^ + i 

z ' 2 


which approaches 1 as z' approaches a pure real, x' < R, since the 
first term in the bracket approaches a pure imaginary. For n odd. 


R.P. 


il 

^sin“- L z > 


d9 = R.P. (9) 


sin' 


-I-* 
z ' 


which may be evaluated as follows: 


Let 


9 = p + iq. 
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then 

sin 0 = sin(p + iq) = sin p cosh q + i cos p sinh q = — 

z ' 

•D 

Since — approaches a pure real as z' — *x', the imaginary term must 
z 1 

be zero; that is, cos p = 0, or p ■ rt/2. Accordingly, 


R.P.(e) 


I sin 


0 


-1-S- 


rt /2 


as z' approaches a pure real, x' < R. The final result for the down- 
wash velocity at the real point z' = x' is, for n odd. 


Ax' n V 

n(n - 2) (n - 4 ) . . . (l) 

2 

Jn - l)(n - 3) . . .(2)_ 


and, for n even, 


Ax ' n V 

n(n - 2) (n - b ) . . . (2) 

2 

{n - l)(n - 3) . . . (1) 


AV Ax * TfV 

For n = 0, the result is for n = 1, it is — jj- — ; for n = 2, it 

is VAx'^j for n = 3, it is ^rtVAx'^. 

8 

Conclusions for Part II 

A study of the symmetry relations concerning the induced velocities 
in the plane of a uniformly loaded rotor indicates that: 

1. The sum of the induced velocities at points (x,y) and (-x,y) 
in the rotor disk is constant and equal to twice the induced velocity 
at the center. 

2. The induced velocity is constant all along the lateral diameter 
of the disk. 

3. The sum of the induced velocities at points (x,y) and (-x,y) in 
the plane of, but outside, the rotor disk equals the induced velocity at 
the corresponding point near the far wake. This sum can be determined 
from the two-dimensional flow about the ellipse that represents the wake 
cross section. 
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A study of the relation between the radial load distribution and the 
radial downwash distribution in the wake of a nonuniformly loaded rotor 
shows that, if the disk loading varies as the nth power of the radius, the 
induced velocity in the far wake also varies as the nth power of the 
radius . 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va., February 24, 1956. 
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Typical load 



Figure 3 *- Disk loadings considered in calculations 
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Figure 4.- Lines of constant value of induced velocity ratio v/vq in 
the longitudinal plane of symmetry of a rotor with a triangular disk 
load. 
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Figure k.- Continued. 
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Figure 4.- Continued. 
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Figure k.- Concluded. 
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Figure 5*- Lines of constant value of induced velocity ratio v/vq in 
the longitudinal plane of symmetry of a rotor with a typical measured 
mean disk load. 
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Figure 5 .- Continued. 
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Figure 5» - Concluded 
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(c) X = ^-5° = tan"- 1 ' 1. 

Figure 6.- Induced velocity distribution along the X-axis of a rotor 

with a triangular disk load. 
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(d) X = 63 . 43 ° = tan -1 2. 



(e) X = 75 . 97° = tan" 1 4. 



(f) X = 90° = tan - -* - 00 . 


Figure 6.- Concluded. 
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Figure 8.- Induced velocity distribution on the Y-axis of a rotor with 

a triangular disk load. 
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Figure 11.- Effect of disk load distribution on the distribution of 
induced velocity along X-axis. X = 63 .^ 3 ° = tan"-*- 2. 
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Figure 12.- Effect of disk load distribution on the distribution of 
induced velocity along the Y-axis. X = 90° = tan“l 00 . 







Figure 13.- Symmetrically located points in the rotor disk. 




Figure 15 . - Two-dimensional crossflow around an ellipse. 
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NACA - Langley Field, Va. 


Figure l6.- Two- dimensional flow around a straight line 


